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INFLUENCE OF MAGNON-PHONON COUPLING ON THE LOW-TEMPERATURE 
MAGNETIC PROPERTIES OF AN ANTIFERROMAGNET* 


P. Pincus and J. WinterT 


Department of Physics, University of California, Berkeley, California 
(Received September 8, 1961) 


There have been several experiments** bearing 
on the temperature dependence of the sublattice 
magnetization, parallel susceptibility, and nucle- 
ar relaxation rates of antiferromagnets. The sur- 
prising result is that, with the possible exception 
of MnF,,° none of the experiments show any evi- 
dence for the expected energy gap® in the spin-wave 
spectrum. Furthermore, the experiments do not 
agree with the predictions of the usual spin-wave 
theory®’” even when the anisotropy is negligible 
and no energy gap is expected. The measure- 
ments”’® of the longitudinal nuclear relaxation 
time T, are also in disagreement with the theo- 
retical predictions.* We give here a theory which 
is able to explain most of these results, but in 
terms of one coupling constant which is not too 
well known independently. The central idea is 
that the magnetostrictive terms in the Hamiltoni- 
an produce a magnon component in the thermal 
phonon spectrum, so that the phonons participate 
directly in magnetic processes. 

We define the temperature Tag corresponding 
to the energy gap by kgTjg =My(2H,H,)™, where 
y is the gyromagnetic ratio, H, the exchange 
field, and H, the anisotropy field. For T<T,4r 
the magnetic thermodynamic properties will be 
mainly determined by the admixture of spin-wave 
character into the acoustic phonons. The effect- 
ive density of thermally excited magnons will then 
not exhibit the exponential temperature variation 
predicted by ordinary spin-wave theory.° 

For the bilinear coupling between the magnons 
and phonons we take 

He’ =G)5.(S*5 76 ** +55 7%), (1) 
Toe =. Se ze 


where G is the coupling constant and €;4” the 


strain tensor. In a forthcoming paper® the trans- 
formations to the mixed magnon-phonon modes 
are given in detail. We have calculated the mag- 
netic properties arising from the low-lying quasi- 
phonons for T<T,g. The deviation of the mag- 
netization is 

2) 2/3 4202 
[a (0) -M(T)]/M (0) = sora a (T/O)'T (484), 

A B (2) 


where M is the ionic mass; a is the lattice con- 
stant; I(x) and ¢(x) are, respectively, the gam- 
ma function and Riemann zeta function; and @ 
is the Debye temperature. For ® =110°K, a=10~ 
em, H,4 =10° oe, and M =10~ g, (2) agrees with 
the T* dependence of the sublattice magnetization 
observed by Poulis and Hardeman’ in CuCl,-2H,O 
if we set G=107** erg. In MnF,, Jaccarino and 
Walker® find that the observed sublattice magne- 
tization does not seem to obey a unique power 
law, but (2) gives the observed order of magni- 
tude when G=107* erg. In MnF,, HA=10* oe, 
@=450°K, and a=3x10~ cm; the other param- 
eters are as for CuCl,-2H,O. This contribution 
to the magnetization deviation is comparable with 
the contribution arising from the ordinary mag- 
non branches of the spectrum in MnF,. 

For T<Ty, x the parallel susceptibility is pro- 
portional to the deviation of the sublattice mag- 
netization and is given by 


xy = (4yh/Ha°)[M(0) -M(T)Y/M(0). (8) 





The parallel susceptibility also has a T* temper- 
ature dependence. Such a variation has been ob- 
served by Stout and Matarrese* for T«T 4g in 
CoF, and FeF,; here G=10~™ erg fits the data 
quantitatively. 
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The low-lying quasi-phonon branch of the spec- 
trum also plays an important role for the longi- 
tudinal nuclear relaxation time T,. We calculate 
the direct relaxation rate, 


~ 31(617)*RAtS*G? | 


1/T, — Ma, OF iy PF /°) (4) 





where A is the hyperfine coupling constant. In 
MnF,, where A =10~ erg, a relaxation rate 
proportional to T has been observed® and agrees 
in order of magnitude with (4) when the param- 
eters above are used. Thus the nuclear relaxa- 

_ tion and magnetization data find a consistent ex- 
planation. The Raman relaxation rate has a T’ 
temperature dependence and agrees again with 
the observations? in CuCl,-2H,O where G=107*° 
erg. The expression for the Raman relaxation 
rate is 





~ 27°(607)¥°A7S°G* | - 
1/T, ro 16RyH ,)"M"a%(k est /®) P,; (5) 


where A is the strength of the dipolar coupling 
and P, is given by 


* max xoo* 
= a a =h / ; 
p,=f = dx; — vk ax’ "pt 
9 (e -1) 


(6) 


Thus by the use of the magnon-phonon interac- 
tion we are able to explain several experiments, ' 
We require, however, somewhat larger coupling 
constants than have been reported in more sym- 
metrical crystals. 

The authors would like to express their appre- 
ciation to Professor C. Kittel for several valua- 
ble discussions. 
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RELATION BETWEEN THERMAL CONDUCTIVITY DEFECT IN SODIUM CHLORIDE 
AND ACOUSTIC RELAXATION EFFECTS AT LOW TEMPERATURES 


K. Brugger and W. P. Mason 


Bell Telephone Laboratories, Murray Hill, New Jersey 
(Received August 3, 1961) 


Recently an anomalously large reduction in 
thermal conductivity at low temperatures has 
been measured’ for NaCl single crystals grown 
by the Harshaw process. No such reduction oc- 
curs for crystals grown free from atmospheric 
contamination, and evidence is presented that 
the reduction in thermal conductivity is associ- 
ated with an impurity ion in the crystal in which 
oxygen is one of the components. The thermal 
conductivity varies nearly as the square of the 
temperature from 1°K to 20°K. Evidence is pre- 
sented that the mean free path satisfies an equa- 
tion of the type 


A =(3.3 107% w/V)* cm, (1) 


where V is the average sound velocity in the De- 
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bye sense and w =27 times the frequency of the 
phonon. Several mechanisms for this effect were 
considered but no conclusive evidence was found 
for any of them. 

One possibility was that a relaxation process 
might be involved. Recently a mechanical relaxa- 
tion effect has been measured at low temperatures 
which is associated with the same set of impurity 
ions, and it is the purpose of this note to show 
that these relaxations are in agreement with the 
loss in thermal conductivity. 

The crystals were obtained from Harshaw and 
from Isomet and were provided with oriented, par 
allel, flat faces. X-cut and AC-cut quartz plates 
were bonded to the specimens with 4-methyl-1- 
pentene. The temperature was controlled by 
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FIG. 1. Attenuation in Harshaw sodium chloride 


crystal for 57-Mc/sec shear wave along [110] with 
polarization direction [110]. Short-dash line shows 
assumed background loss (due mostly to dislocations). 
Long-dash lines show three relaxation mechanisms. 


means of a liquid helium bath and an electric 
heater, and was determined from vapor pressure 
and carbon resistor readings. The attenuation 

was measured with a pulse-echo method. The 
relative velocity was determined using pulse over- 
lap and a phase-comparison technique. 

Figure 1, solid line, shows the attenuation in 
nepers per cm along the [110] direction for a 57- 
Mc/sec shear wave with its polarization in the 
[110] direction. A large increase is seen below 
5°K. No such increase was found for Isomet crys- 
tals, which are grown in a vacuum, as contrasted 
with Harshaw crystals, which are grown in air. 

It is, therefore, concluded that the observed re- 
laxation effect is caused by a gaseous impurity 
introduced from air into the melt.” Similar con- 
clusions on thermal conductivity were found by 
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FIG. 2. Temperature shift of attenuation peak with 
frequency. Inf=Infy-E/kT. 


Klein.’ Table I gives the number of {110} planes 
in which shear stresses are set up by the types 
of waves employed in this investigation and indi- 
cates the corresponding presence or absence of 
a peak. The correlation is excellent. Incidentally, 
the {110} planes are also the slip planes for dislo- 
cations. But in view of the small activation ener- 
gy, a Bordoni relaxation effect can be discounted. 
The temperature shift of the peak as a function 
of the frequency was measured over the range 
from 18 to 173 Mc/sec using shear waves polar - 
ized along [110] and propagating along the [110] 
direction. The result, plotted in Fig. 2, satisfies 
the Arrhenius relation, 


f=f,exp(-E/RT), (2) 


which is typical of a relaxation process. f is the 


Table I. Wave types, shear stress, and attenuation in Harshaw NaCl. 














Number of 

Wave Direction of Direction of {110} planes Attenuation 

mode propagation polarization with shear stress peak 
Longitudinal [100] 4 Yes 
Transverse [100] in (100) 0 No 
Longitudinal [110] 4 Yes 
Transverse [110] {1T0) 6 Yes 

[110] (001) 0 





Transverse 
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measuring frequency and T the peak temperature. 
The constants f, and the activation energy E are 
graphically determined to be 


fo =5X10%*°? cps, E=21+2.3 cal/mole. (3) 


For the same shear wave, the dispersion of the 
sound velocity was measured at 57 Mc/sec ina 
Harshaw crystal as a function of the temperature. 
The curve obtained is typical of a broad relaxa- 
tion process (at 1.5°K the value of V/V, is 0.9986, 
where V is the measured velocity and V, is the 


velocity at the broad flat maximum which stretches 


from 15 to 25°K).° 
A similar attenuation curve was found for KCl 
at 56 Mc/sec with one maximum occurring at 5°K 
and a rise occurring down to 1.5°K. Unfortunately 
no thermal conductivity data are available for this 
crystal. The attenuation peaks of both crystals 
show some evidence of a fine structure and the 
width of the peaks and of the dispersion region 
indicate the existence of a distribution of activa- 
tion energies, rather than a single one. As shown 
by the highest temperature long-dash relaxation 
curve of Fig. 1, a single relaxation agrees with 
the measurement above the peak temperature if 
we assume a background loss— mostly due to dis- 
locations—as shown by the short-dashed line. 
However, at low temperatures, the indicated 
attenuation requires a distribution of relaxation 
frequencies. While the distribution may be a con- 
tinuous one, a good approximation for a slowly 
varying attenuation is obtained by a series of re- 
laxation energies which vary by a small ratio or 
increment.* The attempt frequency of 5x10° cps 
indicates a low probability of occurrence for the 
effect. A relation of the form 
12 (AS;/R) , -(H;/RT) 


f,=10 
i 


(4) 


is assumed for the relaxation frequency of each 
process. The value of 10"? cps is a commonly ob- 
served value for relaxation processes. For the 
highest energy value, AS has to be negative and 
equal to 10.6 cal/mole. The lower the activation 
energy H;, the lower the disorder so that AS; is 
assumed proportional to H;. Then Eq. (4) becomes 


12, -0.252H;,-Hi/RT 


f,=10 (5) 


The internal friction of a series of relaxation 
processes, and the shear modulus change associ- 
ated with them, can be written in the form 


A(f/f.) Au 
Q*s eas — 
1+ (4/7, mM 


i 
loon 
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where A; is a series of constants required to fit 
the shape of the attenuation curve which is related 
to the internal friction curve by the equation 


Qtw if ALf/f) 
a= =>" (7 
2vV V 1+(f/f, 


By assuming activation energy values varying by 
ratios of 1.23,‘ a series of constants A; were ob- 
tained for fitting the attenuation curve of Fig. 1 
above 0.1°K. These constants were tested on the 
available data up to higher frequencies and re- 
produced the measured data with good precision. 
The three top relaxation mechanisms are shown 
by the long-dash curves of Fig. 1. 

Thermal conduction occurs by means of quan- 
tized longitudinal and shear waves in the crystal 
of much higher frequency and these can interact 
with the relaxation mechanisms found by the ultra- 
sonic measurements. Since these phonons come 
from all directions, on the average they will have 
a much smaller component of shearing strain in 
the {110} planes than do the waves whose attenua- 
tion is shown by Fig. 1. A calculation indicates 
that the factor is § for shear phonons which carry 
most of the energy. While a complete calculation 
would require an integration of the distribution of 
frequencies occurring at each temperature, a 
fairly good indication is obtained by taking the 
dominant frequency of the distribution given by 


hf =kT or f =2x10°°T, (8) 


where h is Planck’s constant and k is Boltzmann’s 
constant. Hence, one can calculate the average 
distance required to cause the phonon intensity 
to be reduced to 1/e of its original value —which 
is one definition of the mean free path—by calcu- 
lating from Eq. (7) the distance required to pro- 
duce an attenuation of one neper at the frequency 
corresponding to the temperature T in Eq. (9). 
Using the measured values of a; shown by the 
short-dashed curve of Fig. 1, and multiplying by 
the directivity factor 3, the mean free path is 
shown by the solid line labelled 1 of Fig. 3. The 
shape of the attenuation curve at low temperatures 
is not very critical and the solid line of Fig. 3 
labelled 2 indicates the result of a calculation for 
an attenuation which goes from zero at 0.1°K to 
the measured value at 1.4°K along the dot-dash 
line of Fig. 1. Either curve gives good agreement 
with the thermally measured value shown by the 
dashed line of Fig. 3, since there is no assurance 
that the same number of impurities occurs in each 
sample. The slope of the solid lines varies nearly 
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FIG. 3. Mean free path of thermal phonons. Solid 
line 1 calculated from dashed attenuation curve in Fig. 
1. Solid line 2 calculated from dot-dash curve in Fig. 
1. Dashed line represents Eq. (1). 


as T~'-?8 which is actually in better agreement with 
the measurements? than the assumed slope of T™. 


Thermal resistance caused by a relaxation ef- 
fect appears to be a new process not covered by 
any of the ordinary mechanisms.* It appears that 
the change in the equilibrium states associated 
with the acoustic relaxation effect, produces con- 
siderable attenuation for thermal waves and a 
corresponding short mean free path. Since the 
cross section of the scattering mechanism can 
be obtained from the mean free path and the den- 
sity of scatterers,’ it appears that the scattering 
cross section for relaxation processes of this type 
is much larger than the size of the impurity. The 
very low activation energy range suggests a spec- 
trum of normal modes of low-energy separation 
such as might occur for a tunneling process. 

The authors wish to thank L. M. Tornillo for his 
help in preparing the specimens, Dr. W. Kanzig 
for discussions of the nature of impurities in so- 
dium chloride, and Dr. T. Kjeldaas and Dr. T. 
Geballe for discussions on the characteristics 
of thermal propagation. 





‘Miles V. Klein, Phys. Rev. 122, 1393 (1961). 

{Ww Kanzig and M. H. Cohen, Phys. Rev. Letters 3, 
509 (1959), find O,” ions in NaCl treated with oxygen. 
W. K4nzig (private communication) has evidence for 
the presence of water in the NaCl lattice. N. R. Whel- 
len and T. A. Vanderslice, Bull. Am. Phys. Soc. 6, 
360 (1961), found CO, CO,, and H,O to be given off 
when cleaving alkali halide crystals. 

3In Isomet crystals the sound velocity below 25°K is 
constant to within 1 part in 10°. 

4See Warren P. Mason, Physical Acoustics and the 
Properties of Solids (D. Van Nostrand Company, Prince- 
ton, New Jersey, 1958), Chap. X, p. 300. 

5p. G. Klemens, in Solid State Physics, edited by F. 
Seitz and D. Turnbull (Academic Press, Inc., New 
York, 1958), Vol. 7, pp. 1-98. 
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ENHANCEMENT OF SUPERCONDUCTIVITY BY EXTRACTION OF NORMAL CARRIERS 


R. H. Parmenter 
RCA Laboratories, Princeton, New Jersey 
(Received August 1, 1961) 


One can ask the question: Why does a sufficient- 
ly high temperature always lead to quenching of 
superconductivity in a metal? An examination of 
the Bardeen-Cooper -Schrieffer (BCS) theory of 
superconductivity' shows that the answer to this 
question is as follows: The phonon-induced elec- 
tron-electron attraction between electrons leads 
to a resonant, or coherent, interaction between 
the one-electron states of the metal. The greater 
the number of states interacting, the greater will 
be the resonant interaction, and thus the greater 
the lowering of the energy of the assembly of 
superconducting electrons of the metal. As the 
temperature is raised, electrons are thermally 
excited out of this sea of superconducting elec- 
trons, thus giving rise to normal electrons and 
holes. The one-electron states occupied by these 
normal carriers are no longer available for the 
resonant interaction (because of the exclusion 
principle), thereby decreasing the strength of the 
resonant depression of the energy of the system. 
At a sufficiently high temperature (the supercon- 
ducting transition temperature), there are suf- 
ficient normal carriers present to quench the 
resonance completely, whereupon the metal be- 
comes normal. 

This suggests that, if there were some method 
of electrically extracting the normal carriers as 
soon as they are thermally generated, then it 
might be possible to maintain the metal super- 
conducting above the usual transition temperature. 
As a result of recent work,’ it is known that it is 
indeed possible to inject or extract normal car- 
riers from a superconductor. (It is this phenom- 
enon which insures the success of the operation 
of the superconducting tunnel diode.) This is 
done by making contact to another superconductor 
through a thin (10-20 angstroms) insulating layer, 
the latter serving to pass normal carriers under 
suitable bias (by tunneling), but to block super- 
conducting electrons completely. It appears that 
an insulating barrier more than a few angstroms 
thick is effective in preventing current flow in- 
volving the cooperative motion of many electrons 
(i.e., the superconducting current) while at the 
same time allowing some current flow involving 
uncorrelated motions of individual carriers (i.e., 
the normal current). 

Let us suppose we wish to extract normal car- 
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riers from superconductor A having transition 
temperature 7,4. We assume that the operating 
temperature T is close to (but smaller than) T oA 
in order that there be appreciable numbers of 
normal carriers in A under equilibrium condi- 
tions. Through a thin insulating layer, we make 
contact with superconductor B having transition 
temperature 7,p sufficiently greater than 7 for 
there to be negligible numbers of normal carriers 
present in it. By now biasing such a diode with a 
voltage equal to one-half the energy-gap difference 
of the two superconductors, we can extract either 
normal electrons or normal hoies (depending on 
the polarity of the bias) from A. By having A in 
the form of a film sandwiched between contacts 
(electrodes) of superconductor B (along with the 
thin insulating layers), we can extract both signs 
of normal carriers by putting a bias between the 
two B contacts equal to or greater than the energy- 
gap difference 2(egp-€ 9,4). If the bias is made 
greater than the energy-gap sum 2(€gp+ €g,), 
however, we get an additional effect due to gener- 
ation of isoenergetic electron-hole pairs at each 
insulating interface—this effect leading to a double 
injection of normal carriers into A. This injec- 
tion apparently will overbalance the previously 
discussed double extraction, so that for biases 
greater than 2(€gp+€Q,), there will be a net in- 
jection rather than extraction. Figure 1 illus- 
trates the energy-level diagrams for such a device 
under conditions of zero, 2(€gp-€g,4), and 
2(€9R+ €QA4) bias. 

We designate by 7 the lifetime against normal 
electron-hole recombination and by n’ the densi- 
ty of normal electrons (or holes) in a supercon- 
ductor under thermal equilibrium conditions. 
Thus n’/t is the normal electron-hole pair gen- 
eration rate per unit volume in the superconduc- 
tor. (It is also the recombination rate.) Under 
conditions of extraction in the superconductor, we 
shall, for the moment, assume the normal pair 
generation rate is unchanged (i.e., that it remains 
n'/t). The recombination rate, however, will be 
reduced to (n”/n’)?(n’/1), where n” is the reduced 
density of normal electrons resulting from the ex- 
traction process. The factor (n”/n’)? is a conse- 
quence of the fact that we are dealing with a bi- 
molecular recombination process. We designate 
by L the thickness of the superconducting film 
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FIG. 1. Energy level diagrams for a system of a 
superconducting film of thickness L and energy gap 
2€9,4 Separated by insulating layers from superconduct- 
ing contacts having energy gap 2€gg. (a) Zero bias. 

(b) Bias = 2(€ OB ~ € 0A) (c) Bias = 2(€ opt € 0A)* 
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being subjected to extraction, by y the mean tun- 
neling transition probability per collision of the 
normal carriers striking the insulating barriers, 
and by vy the velocity of the normal carriers 
(this being the Fermi velocity). The pair extrac- 
tion rate per unit volume of superconducting film 
is jup(n”/2)(y/L). Upon setting the difference 
between the generation rate and the recombina- 
tion rate equal to the extraction rate (i.e., dy- 
namic equilibrium), we get 


-G)-@EIE) © 


where \=vp-7 is the mean free path against pair 
recombination under thermal equilibrium condi- 
tions. If yA/6L>1, then 


n"/n'=6L/yA«1, (2) 


and the extraction process is very efficient. An 


order -of-magnitude lower limit to A is given by 


the normal-state bulk electrical conductivity mean 


free path A, as limited by lattice vibrations, d,, 
being a lower limit since it is harder to generate 
normal pairs thermally in a superconductor than 
in the corresponding normal metal. The real 
phonon absorbed in the generation process must 
have at least the gap energy in the case of the 
superconductor. Thus if the operating tempera- 
ture is appreciably below the usual transition 
temperature of the superconductor, we can ex- 
pect A>>A,.- Since r,, may be 10™ cm and y 
~1077 - 107° under typical conditions, Eq. (2) 
shows that L, the superconducting film thickness, 
may be 100 angstroms or more and still have 
(n”/n')«1. 

In carrying out the analysis of the extraction 
process, we have assumed spatial uniformity of 
n” in the superconducting film. This will be an 
excellent approximation as long as L«). We 
also assumed the pair generation rate to be un- 
affected by the extraction process. As we shall 
see in a moment, however, the superconducting 
energy gap may increase under conditions of ex- 
traction. If so, this implies a decrease of the 
generation rate with increasing extraction (fewer 
phonons are energetic enough to create pairs). 
Such a decrease of the generation rate obviously 
does not invalidate the conclusions of the pre- 
ceding paragraph. 

A third assumption is that the density of nor- 
mal carriers in the contacts is negligible com- 
pared with n”, the density in the superconducting 
film, so that there is negligible tunneling of nor- 
mal carriers from the contacts back into the film. 
This implies that there is an efficient method of 
getting rid of the normal carriers injected into 
the contacts. One way of accomplishing this is to 
plate each contact with a certain amount of nor- 
mal metal. The latter will serve as an efficient 
sink for normal carriers injected into the con- 
tact, whenever the distance from the tunnelable 
barrier facing the contact is less than a mean 
free path (for normal carriers in Sz) from the 
normal metal of the contact. The normal metal 
should be separated from the tunnelable barrier 
by a thickness of superconductor B of at least 
something of the order of the Pippard coherence 
distance (10 cm), and the normal metal itself 
should have minimum thickness of similar size. 
This obviates the possibility that the normal 
metal is made superconducting by being in con- 
tact with Sp, or conversely that Sp has its super- 
conductivity quenched by being in contact with a 
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normal metal.® 

It should be pointed out that the superconductor 
Sp in each contact contains injected normal elec- 
trons or normal holes, but not both. (One con- 
tact will have electrons, and the other holes.) 
The resulting space charge in each contact will 
be compensated by an adjustment in the density 
of superconducting electrons (a shift of the ener- 
gy of the bottom of the conduction band in S$ B rela- 
tive to the Fermi level). This will lead to a mi- 
nute shift in the transition temperature of each 
contact.* Such a shift will be negligible compared 
with the shift of transition temperature of the 
superconducting film S, due to extraction of nor- 
mal carriers of both signs without change of the 
density of superconducting electrons. 

The situation we are considering in this note 
represents a nonequilibrium condition. There 
are two general types of nonequilibrium condi- 
tions in a superconductor: (1) those where the 
lack of equilibrium is due to the superconducting 
electrons, so that dissipation processes do not 
occur; (2) those where the lack of equilibrium is 
due to the normal carriers, so that dissipation 
processes do occur (i.e., finite entropy produc- 
tion). Time-independent current flow in a super- 
conducting wire is an example of the first type.°® 
Because of the lack of dissipation, it is possible 
to define a free energy, despite the lack of equili- 
brium. Specifically, one modifies the equilibrium 
free energy by adding a term equal to the product 
of the quantity being constrained times a Lagrang- 
ian multiplier. (The constraint, set by boundary 
conditions, is what causes disequilibrium.) For 
the example of dc current flow in a superconduct- 
ing wire, the constraint is that imposed on the 
net current due to the superconducting electrons.® 

The extraction of normal carriers represents 
a nonequilibrium situation of the second type 
where dissipation processes occur. There is a 
certain amount of thermoelectric cooling in the 
superconducting film being exhausted of normal 
carriers, and a still greater amount of heating 
in the contacts. For example, under the bias con- 
ditions of Fig. 1(b), heat is being removed from 
the film at the rate J(2€9,/e) and being liberated 
in each contact at the rate I(€gp/e) when extrac- 
tion is efficient. (Here J is the total current, and 
e is the electronic charge.) The removal of heat 
in the film occurs when phonons are absorbed in 
making normal hole-electron pairs. Because of 
backflow of heat from the contacts into the film, 
the actual temperature drop of the film is probab- 
ly small enough to ignore. Nevertheless, the 
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presence of dissipation associated with this ther- 
moelectric process makes it diticult to define 
any sort of free energy F such that the minimiza- 
tion of F will lead to a description of steady- 
state equilibrium. 

A generalization of the BCS theory to such a 
nonequilibrium situation can be made in the follow- 
ing manner. A Boltzmann transport equation is 
set up for the distribution function f, associated 
with the normal carriers (normal electrons for 
k>kp, and normal holes for k<kp, k being the 
wave vector labeling the single-particle states, 
and k ;, being the Fermi wave vector). At the 
same time, the internal energy U is minimized 
with respect to the parameters hp appearing in 
the BCS many-electron wave functions.” Such a 
minimization leads to the equation 


hw 
13. » 

[v(0)V] {f (c,2+e 277 (1-2, Mde,. — ) 
Here N(0)V and fw are constraints characteristic 
of the superconductor, €, is the single-electron 
Bloch energy (measured with respect to the Fermi 
level), 2€, is the energy gap, and (€,?+€ 9") is 
the energy of a normal carrier in the supercon- 
ductor (also measured with respect to the Fermi 
level). Equation (3) is formally the same as that 
of the BCS theory; the only difference lies in the 
fact that the f, contained therein is obtained from 
a Boltzmann transport equation rather than from 
a minimization of a free energy with respect to 
Sr: 

Of course, it is possible to have a nonequili- 
brium situation where both the superconducting 
electrons and the normal carriers contribute to 
the disequilibrium. For such a situation, one 
should solve a Boltzmann equation for fp and at 
the same time minimize the internal energy U 
with respect to hp subject to the constraint on the 
superconducting electrons causing their disequili- 
brium. Such a situation arises in the particular 
physical problem of the superconducting film with 
extracting electrodes if we choose to make a su- 
percurrent flow in the film parallel to the plane 
of the film at the same time that we are extract- 
ing normal carriers from the film. 

For the moment, we return to the case where 
there is no supercurrent in the superconducting 
film. Boltzmann’s equation for our problem re- 
duces to the statement that the total time rate of 
change of f, vanishes, the three processes of 
normal-pair thermal generation, normal-pair 
recombination, and double extraction of normal 
carriers all contributing to (df,/dt). The extrac- 
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tion process is effective in lowering f, for all 
orientations of k, despite the fact that only car- 
riers moving nearly normal to the insulating in- 
terface have appreciable probability of tunneling. 
This lowering of f, for all orientations is a con- 
sequence of the fact that those carriers reflected 
at the insulating interface are undoubtedly re- 
flected diffusely, thus leading to rapid equilibra- 
tion between isoenergetic carriers moving in 
different directions. 

It is not really necessary to solve Boltzmann’s 
equation for f,. Under reasonable physical con- 
ditions, it appears that whenever extraction is 
efficient as a whole, i.e., whenever (n” /n’)«1, 
it will also happen that 


hi. (4) 


This will hold for all k if the bias is adjusted such 
that normal carriers of all energies can tunnel 
[e.g., as in Fig. 1(b)]. [In contrast, under no- 
extraction conditions where n’ is appreciable, 
there will be a number of f, comparable with 5 

and thus not satisfying (4).] It should be noted 

that the only way the temperature T can affect 

the energy gap 2€,, obtained by solving Eq. (3), 

is through the temperature dependence of f,. How- 
ever, as long as (4) holds true, then (3) will be 
independent of T, and we will obtain an €, equal 

to the BCS value at the absolute zero of tempera- 
ture, even when T is above the usual supercon- 
ducting transition temperature. 

It would appear that, by means of electrical ex- 
traction, we can raise the transition temperature 
of the superconducting film. An upper limit to 
this new transition temperature will, of course, 
be set by T, p, the transition temperature of the 
contacts. As the operating temperature approaches 
Tp, the extraction efficiency will drop off be- 
cause of the presence of thermally generated nor- 
mal carriers in the contacts which may tunnel 
through the insulating layers into the supercon- 
ducting film. It appears plausible that one should 
be able to start with the film in the normal state 





and by means of extraction cause the film to go 
superconducting, the operating temperature being 
less than T, p but greater than T,, (the usual 
transition temperature of the film). Conversely, 
we should be able to lower the transition tempera- 
ture of the film by injection of normal carriers 

as shown in Fig. 1(c). By lowering the transition 
temperature below the operating temperature, we 
can drive the superconducting film normal by elec- 
trical injection. Using either injection or extrac- 
tion, we have an electrical (nonmagnetic) method 
of controlling a superconducting current in the 
film (the current flowing parallel to the plane of 
the film). 
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®In reference 5, it appears that the constraint is being 
imposed on the normal current rather than the super- 
conducting current. This is because the calculation, 
for convenience, is being carried out in a coordinate 
system moving with the mean drift velocity of the super- 
conducting electrons. In the laboratory coordinate sys- 
tem, the constraint is being imposed on the supercur- 
rent, not the normal current. 

'The internal energy U is formally the same as in the 
BCS theory, i.e., Eq. (3.16) of BCS after removal of 
the term -7S. In the BCS theory, the free energy 
F=U-TS, rather than U, is minimized with respect to 
h,. This difference is of no consequence, since the 
entropy S is independent of hy. The minimization of 
U with respect to hp, appears to be the only reasonable 
way of generalizing the BCS variational approach to a 
nonequilibrium situation. Such a minimization never- 
theless represents an assumption. 
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SUPERCONDUCTIVITY OF SUPERIMPOSED METALS* 


A. C. Rose-Innesf and B. Serin 
Rutgers University, New Brunswick, New Jersey 
(Received September 1, 1961) 


Measurements? * on thin films of superconduc- 
tors in contact with normal metals have shown 
that the transition temperature of the films is low- 
ered. For sufficiently thin films the transition 
temperature may be reduced to such an extent that 
superconductivity is not observed. There have 
been several theories of this effect.'»*»° Our ex- 
periments described here indicate that for films 
deposited by evaporation the effects observed can 
be considerably influenced by diffusion between 
the two metals. 

We have made a series of measurements of the 
magnetic moment of superconducting tin films 
evaporated onto gold plates. The method used to 
detect the magnetic moment was similar to Alek- 
seyevsky’s,® in which the specimen is suspended 
from a torsion balance at a small angle to a uni- 
form magnetic field. As the field strength is in- 
creased the film rotates until it is parallel to the 
field. The deflection is then constant until the 
critical field value is reached, when the specimen 
returns to its original position. This method is 
preferable to the measurement of resistance which 
has been used in all previous experiments. The 
resistance method has the disadvantage that it 
may not necessarily measure a bulk property of 
the film. Zero resistance may arise from super- 
conducting filaments, or a finite resistance may 
result from small breaks in the film, neither of 
which is characteristic of the material as a whole. 
However, the magnetization is a bulk property, 
not affected by small discontinuities. Moreover, 
measurements of the magnetic moment as a func- 
tion of temperature and magnetic field show wheth- 
er or not the specimens behave like films of an 
ideal superconductor. 

The substrates for our films were 99.99% pure 
gold plates 21 mm x 21 mm and of thickness 0.25 
mm. The plates were annealed at 750°C for about 
a day and then wiped with trichlormethylene. Dur- 
ing the evaporation they were firmly clamped to a 
hollow copper block cooled by flowing water at 
room temperature. The tin was evaporated from 
a molybdenum boat about 15cm below the plate. 
During evaporation the pressure was less than 
5x107° mm Hg. At the same time, the film was 
also deposited onto a thin glass cover slip which 
was next to the gold plate and symmetrically 
placed with respect to the boat. Between the boat 
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and the plates was a shutter which was opened only 
after the tin had been melted, outgassed, and 
evaporation had begun. The shutter was kept open 
for 2.5 min for all evaporations so that every 
specimen saw the hot boat for the same length of 
time. The electrical power supplied to the boat 
also was the same for all depositions. 

The anount of tin placed in the boat was weighed 
to give the thickness desired, and the evaporation 
time was sufficient to insure that all the tin vapor- 
ized. Thicknesses quoted in Table I were deter- 
mined by weighing the glass slips before and after 
deposition. The films were squares of area about 
3 cm’. 

Table I shows the results of the measurements 
on these films. There is a rapid decrease in tran- 
sition temperature when the thickness of tin is 
less than about 7000 A; also, the films fall into 
two groups. Those specimens with transition tem- 
peratures close to the value of bulk tin (say T; 
> 3.4°K) had sharply defined critical magnetic 
fields, and the magnetization and critical fields 
exhibited the behavior to be expected for thin films 
of an ideal superconductor. However, those films 
whose transition temperatures were appreciably 
depressed showed the characteristics of nonideal 
superconductors, having very high critical fields 
which were not sharply defined. Furthermore, 
their magnetic moments were smaller than those 
of the other films. The films deposited onto 
glass all had transition temperatures near the 


Table I. Superconducting properties of tin films on 
gold. 








Film Thickness (A) T. (°K) Behavior 
5 18 700 +500 3.71 ideal 
8 13 200 3.40 ideal 
9 9500 3.51 ideal 
10 7800 2.68 nonideal 
13 7100 3.57 ideal 
14 7100 3.59 ideal 
15 6100 3.08 nonideal 
11 6100 1.66 nonideal 
. one not superconducting 
° 7 down to 1.35°K 
4 1500 
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value of bulk tin. 

These results appeared to confirm qualitatively 
the earlier observations of this effect. Misener 
and Wilhelm! found that tin films electroplated on 
constantan ceased to be superconducting above 
2°K at thicknesses less than 3000 A. More recent- 
ly, Meissner? measured several pairs of metals, 
including tin on gold, and reported a critical thick- 
ness of a few hundred angstroms. Smith, Shapiro, 
Miles, and Nicol® observed qualitatively similar 
effects with lead and silver films. 

However, by further experiments we showed 
that under the conditions of our evaporations dif- 
fusion probably occurred despite the cooling pre- 
cautions. This was established by replacing some 
of the specimens in the evaporator and exposing 
the films to the radiation from the hot empty boat, 
under conditions similar to those used in laying 
them down. After this treatment we measured 
their superconducting properties again and found 
them changed. The transition temperature of film 
14 was reduced to 2.26°K after a 2.5-min exposure, 
and its magnetic behavior changed from ideal to 
nonideal in the sense already explained. After a 
further exposure of 2.5 min, film 14 was not su- 
perconducting down to 1.27°K. Film 9 after a 5- 
min exposure was not superconducting down to 
1.23°K. The appearance of the films was unchanged 
after these exposures. Since the superconductivity 
of our films can be destroyed by heating the tin- 
gold surface, it follows that the value of the thick- 
ness at which the transition temperature is ob- 
served to fall rapidly has no particular signifi- 
cance, but is a consequence of the way the speci- 
mens were prepared. 

We believe that the effects of reheating show 
that diffusion can easily occur during evaporation, 
and result in a decreased transition temperature. 
Furthermore, diffusion might account for the dif- 
ferent values given by various workers for the 
thickness at which superconductivity disappears. 
The degree of diffusion and alloying which can 
take place will, of course, depend on the partic- 
ular pair of metals used and the method of film 
deposition. 

Gold is only soluble in tin up to 0.3% and about 
4% of tin is soluble in gold, but intermetallic 
compounds AuSn, and AuSn, exist.’ Allen® studied 


the superconductivity of tin-gold alloys and found 
that the transition temperature of the tin-rich 
alloys is independent of gold concentration up to 
about 30% gold. The intermetallic compounds are 
also superconducting, but alloys with more than 
50% gold are not superconducting down to 1.5°K. 

Although diffusion seems to reduce the transi- 
tion temperature, the exact reason for this is 
not obvious. Appreciable interpenetration of the 
two metals could have two effects: Nonsupercon- 
ducting alloys may be formed; or the thickness 
of the superconducting metal could be reduced to 
a value much less than that deduced from the 
weight. Alternatively, there may be only a small 
interpenetration resulting in more perfect contact 
between the two metals,°® which according to Coo- 
per® can lower the transition temperature of rel- 
atively thick films. 

In conclusion, we feel that we have shown that 
diffusion can play a critical role in experiments 
on the superconductivity of superimposed metals; 
and consequently, great caution must be used in 
the interpretation of such experiments until the 
role which diffusion plays is adequately under - 
stood. 

We are grateful to Miss S. Goldis for her assist- 
ance in these experiments, 
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ABSENCE OF ANTIFERROMAGNETIC DOMAIN WALLS IN MnF, 


P. S. Pershan 
Gordon McKay Laboratory, Harvard University, Cambridge, Massachusetts 
(Received September 5, 1961) 


It has been suggested that the electrically in- 
duced shift of the F*® nuclear resonance frequency 
in MnF, can be used to study the presence or ab- 
sence of antiferromagnetic domain walls.’ The 
purpose of this note is to report the absence of 
domain walls in a single crystal of MnF,. 

In MnF,, below the Néel temperature, the anti- 
ferromagnetic ordering of the Mn ions causes 
a large internal magnetic field at the F sites; at 
4.2°K, with no applied magnetic field, the F nu- 
clei resonate at 159.970 Mc/sec.? 

Figure 1 shows two alternate choices for the 
unit cell of MnF,, depending on which Mn*~ ion 
is chosen for the corner. For the F nuclei at A 
sites, the three nearest Mn‘ ions in (110) planes 
and their magnetic moments are oriented so the 
local field at the F is up. For the B sites the 
three nearest neighbors are in (110) planes and 
the local field is down.® 

An electric field in the [110] direction will split 
the nuclei at A sites by 1.4 kc/sec per volt/cm 
but leave the B sites unaffected.’ A magnetic 
field in the [001] direction, as shown, and no elec- 
tric field, will increase the resonance frequen- 
cies of the nuclei at A sites by 4 kc/sec per gauss 
and decrease the ones at B sites by the same 
amount. If both these fields are simultaneously 
applied, the resonance frequencies of the nuclei 
at A sites will increase and split into a doublet 
while the nuclei at B sites will decrease their 
frequency and remain unsplit. 

If, however, the crystal has another domain, 
the Mn** magnetic moments in that domain will 
be antiparallel to the arrows in Fig. 1. When the 
field H is applied as shown, the resonance fre- 
quency of the nuclei at A sites will decrease rath- 
er than increase. For the electric and magnetic 
fields as shown, this domain would have a doublet 
at a lower frequency than the singlet. Superposi- 
tion of the resonances from both domains would 
result in two triplets; the relative strength of the 
outer to inner components of the triplet is a meas- 
ure of the fraction of crystal taken up by each 
domain. For a multidomain crystal both types of 
domains occupy equal volume and the central com- 
ponent of the triplet is twice as big as the outer 


x7.9 mm with the [110] direction perpendicular to 
the broad face. Electrodes were applied to the 
broad face as described in reference 1. The resgo- 
nances were observed with the spectrometer built 
by Kushida for the pressure dependeace of the F 
resonance in MnF,.* The spectrometer was fre- 
quency swept by a motor-driven capacitor and 
frequency modulation was used to record the first 
derivative of the lines. It was not possible to oper- 
ate the spectrometer without some degree of sat- 
uration; hence one must exercise caution in try- 
ing to interpret linewidth and amplitude variations 
of the traces in Fig. 2. The magnetic field was 
supplied by a small permanent magnet. 


The resonance in Fig. 2(a) is with no external 
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FIG. 1. Two alternate choices for the unit cell of 








MnF,, depending on the particular Mn** ion chosen 
for a corner position. The arrows indicate the orien- 
tation of the time-averaged magnetic moment for the 
Mn** ion. 


components. 
The crystal was cut from a single crystal of 
MnF, as a plane parallel slab 1 mmx4.8 mm 
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FIG. 2. The F’* nuclear resonance signal with vari- 
ous combinations of electric and magnetic fields. 


fields. In Fig. 2(b) the line is split by an electric 
field in the [110] direction; the doublet due to 
spins at A sites combine with a singlet from spins 
at B sites to form atriplet. For Fig. 2(c) the 
magnetic field is on and the electric field is off. 
Figure 2(d) has both electric and magnetic fields 
applied as shown in Fig. 1. If we compare the 
shape of the doublet with the triplet in Fig. 2(b), 
there is no evidence for more than one domain. 
With the signal-to-noise level available, we can 
conclude that more than 80% of the crystal is of 
one domain type. In Fig. 2(e) the magnetic field 
is reversed interchanging the singlet and doublet. 
The singlets in 2(d) and 2(e) seem to be broad- 
ened by the electric field much more than either 
component of the doublets. This is expected 

since the singlets are sensitive to misorientation 
of the electric field to the first order while the 
doublets have only a second-order dependence. 

These same results were obtained on two oc- 
casions. The only difference was that, on warm- 
ing up and recooling below the Néel temperature, 
the two sublattices reversed direction. 

Absence of domain walls in perfect crystals of 
MnF, at 4.2°K is expected on the basis of theo- 
retical calculation of Li.’ The large anisotropy 
forces in MnF, require a domain wall to be only 


one lattice spacing thick. There is no reduction 
in the magnetostatic energy to compensate for this 
large increase in the exchange energy and the gain 
in entropy on creating a domain wall is not large 
enough to make the change in free energy negative. 

Antiferromagnetic domains can, and probably 
do, exist for some time after the crystal is first 
lowered through the Néel point. Well below the 
Néel point they are thermodynamically unstable, 
however, and unless pinned down by some sort 
of lattice inperfection they diffuse themselves out 
of existence. 

The stability and motion of domain walls near 
the Néel point has recently received theoretical 
attention.* Observations of the F*® nuclear reso- 
nance just below the transition temperature show 
linewidths up to 600 kc/sec.” This is much larger 
than the maximum electric field splitting observ- 
able in MnF,. The linewidth, however, is a strong 
function of temperature and it should be possible 
to search for domain walls about 10°K below the 
Néel point. 

This research was stimulated by discussions 
with Professor N. Bloembergen and grateful ac- 
knowledgment is given for his contribution. 
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NATURE OF THE V, CENTER 


Frederick Seitz 
Department of Physics, University of Illinois, Urbana, Illinois 
(Received September 11, 1961) 


The V, center, discovered by Casler, Pring- 
sheim, and Yuster’ in 1950, has been the object 
of much speculation. Operating within the frame- 
work of pure vacancy models, which seemed most 
appropriate a decade ago, the writer? proposed 
that this center is the antimorph of the F center, 
consisting of a positive-ion vacancy and a hole. 

In the meantime, the experimental work of Kan- 
zig, Woodruff, and Castner® and Delbecq, Smaller, 
Yuster, and Hayes* has made it increasingly clear 
that many centers of the V type very probably in- 
volve interstitial halogen atoms, in keeping with 

a general conjecture first put forward by Varley.® 
In particular it seems highly probable that the H 
and V, centers, which are closely related, are of 
this type. Kingsley in connection with his ex- 
tensive study of the centers, has given evidence 

to show that the present author’s model of the V, 
center may well be assigned more appropriately 
to the experimentally observed V, center. 

There is strong support for the view that the H 
center is a neutral halogen atom in a crowdion 
configuration in a (110) row of halogen ions. Pre- 
sumably the center is normally produced in close 
juxtaposition to an F, F’, or a center when gen- 
erated by excitation of the electronic structure of 
the lattice. Kanzig, Castner, and Woodruff*® have 
proposed that the V, center is a neutral halogen 
molecule in a vacant halogen site—a model which 
faces numerous difficulties. Cape’ has proposed 
another model. It seems hopeless at present to 
propose a logically water-tight model of the V, 
center which ties together in a highly precise way 
all of the observations made on different crystals 
of varying perfection in different laboratories by 
different observers. It seems to the author, how- 
ever, that the model which combines the maximum 
degree of simplicity with the minimum number of 
major contradictions is based on the view that the 
V, center is an interstitial negative ion, presum- 
ably surrounded tetrahedrally by like and unlike 
neighbors [i.e., at a (¢,4,4) position]. It is as- 
sumed that this center migrates by interstitialcy 
diffusion when the V, band bleaches so that the Dut- 
ton-Maurer current® is electrolytic. Since the in- 
terstitial ions should be close to halogen-vacancy 
partners initially, it is reasonable to expect® the 
current peak to be displaced from the annealing 
peak for optical absorption, provided the V, 
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centers are annihilated at F or a centers at the 
terminus of migration. It is assumed that the con- 
version from H to V, in warming from liquid heli- 
um to liquid nitrogen range of temperature involves 
the transfer of an electron from an F center to an# 
center to form an energetically more stable vacan- 
cy-interstitial pair. It should be stated that this 
model was mentioned in passing by Kingsley,® but 
was discarded for reasons which possibly are too 
severe. Some of its principal virtues are as fol- 
lows: 

1. The absence of electron spin. 

2. The absence of dichroism.?° 

3. It is negatively charged so that it will not 
capture photoelectrons from F or F’ centers. 

4. F’ centers are immune to annihilation by 
migrating V, centers but F and a centers are not. 
5. The absence of an observable Hall current 

when migration starts. 

6. It is a natural complement to the halogen- 
vacancy centers so that there is no residual inter- 
stitial when it combines with an a or F center. 

7. One may reasonably expect it to migrate 
under repeated optical excitations. 

Apart from questions relating to the relative 
stability of the V, center proposed here, the model 
faces the following major obstacles: 

1, The V, center does not show appreciable 
photoconductivity, as one might naively expect 
from its negative charge and by analogy with the 
F’ center. Actually this appears to be a delicate 
problem since the center has an appreciable 
Coulomb-like polarization field. It is possible 
that the first excited state is stable. 

2. One wonders why the V, center does not 
blunder into some of the V; centers (self-trapped 
holes) when they migrate and annihilate a sig- 
nificant fraction of the latter. If they did, H cen- 
ters would be produced. 

3. Kingsley® points out that if an interstitial 
halogen ion were to meet an F center, a free elec- 
tron which could annihilate a Vz center would be 
produced whereas no such effect appears to be 
observed. This matter remains open for ques- 
tion. 

The model might be tested in KBr by seeing if 
additions of KCl inhibit both thermal bleaching 
and charge flow associated with V,. Presumably 
interstitial Cl” ions, once formed, would trap the 
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interstitialcy pattern. 

There has been frequent discussion" of the 
manner in which excitation of the electronic sys- 
tem may lead to the generation of vacancy-inter- 
stitial pairs in the halogen lattice. It may not be 
out of place to point out that an @ center (a halo- 
gen-ion vacancy) is created when a U center (a 
negative hydrogen ion in the site of a halogen ion) 
is excited with light at low temperatures.” It fol- 
lows that the excited hydrogen ion is sucked into 
an interstitial position as a result of the forces of 
attraction between it and the normal halogen ions 
of the lattice. There seems to be no fundamental 
reason why this should not occur to an excited hal- 
ogen ion under appropriate conditions.’* Thus a 
self-trapped exciton may spontaneously generate 
a halogen-ion vacancy and an interstitial halogen 
ion. In the case of the halogen, the generation of 
the interstitial configuration could involve inter- 
stitialcy motion, the ion initially at the site of the 
vacancy moving to a normal nearest-neighbor po- 
sition. 

The author is indebted to Dr. Liity for a discus- 
sion of this topic. 
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LEWIS EFFECT IN RESONANCE YIELD CURVES* 


W. L. Walters,? D. G. Costello, J. G. Skofronick, D. W. Palmer, W. E. Kane, and R. G. Herb 
Department of Physics, University of Wisconsin, Madison, Wisconsin 
(Received September 5, 1961) 


Lewis’ has recently pointed out that a thick- 
target yield curve taken at a narrow resonance 
with a charged particle beam of high-energy res- 
olution should have a maximum just above the 
resonance energy followed by a shallow minimum, 
before it assumes a constant value. The qualita- 
tive explanation for this phenomenon, which the 
authors have chosen to call “the Lewis effect, ” 
is as follows: In passing through target material 
a charged particle loses energy in discrete steps 
Q. If some of these steps are larger than the 
natural width of a narrow resonance, some of the 
particles incident on a target at an energy well 
above the resonance energy, E R? will jump over 
the resonance. If particles are incident at Ep 
then all will have for a finite time the correct 
energy to interact. The yield curve should there- 
fore exhibit a peak near Ep. 

Gamma-ray yield curves from the Al?"(p, y)Si”® 
resonance reaction at 992 kev were studied with 
a proton beam having an energy spread at half 
maximum of about 125 ev. The natural width of 
this resonance is 80440 ev.” Yield curves were 
initially taken from aluminum films formed in 
the target chamber by evaporation of aluminum 
from a tungsten filament. A peak in the yield 
curves appeared but could not be reproduced. 
This was interpreted as being due to the accu- 
mulation of surface deposits on the target. A 
new method of target preparation was developed*® 
in which aluminum is slowly and continuously de- 
posited on the target backing while taking data. 

Figure 1 shows several yield curves exhibiting 
the Lewis effect. Curve A is a yield curve cal- 
culated as outlined below. The data labelled B 
and C were taken during continuous evaporation 
of aluminum. Curves D and E were initial runs 
from two different targets prepared by filament 
evaporation; subsequent runs in both cases failed 
to reveal the Lewis peak. 

The yield per proton at a mean beam energy E, 
from an infinitely thick homogeneous target is 
given by 


¥é,) “bp aods = 08% 8) 


x o(E ,,E)n(E, E )dEdE ,, (1) 


where n, is the number of aluminum atoms per 
unit volume, g(E,,£;)dE; is the probability that 
a proton in the beam of mean energy Ey» has an 
incident energy between E; and E;+dE;, o(Ep,E) 
is the cross section for the resonance, and n(E, 
E£;)dE represents the probability that a proton 
incident at an energy E; has an energy between £ 
and E +dE somewhere inside the target. 

The probability p(Q, E)dQ that a proton at an 
energy E will suffer a collision in which it loses 
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FIG, 1. Thick-target yield curves showing the Lewis 
effect from the Al?"(p, y)Si®® resonance at 992 kev. Mean 
beam energy, Ep, is plotted relative to the resonance 
energy, Ep. The position of Ep is determined by the 
calculated yield curve. Plateau yields are normalized 
for easy comparison of peak amplitudes. 
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energy between @ and Q+dQ is given by* 
p(Q, E)dQ = (K/EQ”)dQ, (2) 


where K is a constant. Using for this resonance 
Qmin =12-3 ev and Qmax = 2160 ev, the function 
n(E, £;) was numerically calculated on a computer 
using a Monte Carlo method. 

The function n(E, £;) gives the yield per proton 
at energy £; from a resonance of infinitesimal 
width with a monoenergetic beam incident on the 
target. To account for nonzero resonance width 
and finite-beam energy resolution, the products 
of n(£,£;) values and thin-target yield curve val- 
ues were graphically integrated to give curve A 
in Fig. 1. The thin-target curve chosen for the 
calculation was from a target of 60-ev thickness. 

Curve A exhibits the peak and shallow minimum 
predicted by Lewis. The amplitudes of the peaks 
on experimental yield curves vary by a factor of 
two. This strongly suggests that surface contam- 
ination of the target plays an important role even 
on targets being continuously evaporated. Some 
of the discrepancy between the calculated curve 
of Fig. 1 and the data points B, which exhibit the 
highest experimentally observed peak, may be at- 
tributed to contaminants. On the other hand, the 


peak in the calculated curve may also be too high 
because of approximations in stopping-power the- 
ory. 

The point of half plateau yield on the calculated 
yield curve, which is commonly assumed to be 
ER, does not fall at Er. For the parameters 
used here, the half-value point is about 100 ev 
below Ep. 

Experimental work is being extended to other 
resonance reactions and an attempt will soon be 
made to utilize higher beam energy resolution. 
Calculations have shown that the Lewis peak in- 
creases in amplitude as resolution improves and 
as the values of Qmax and Qmin increase. 
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ASYMMETRY PARAMETERS IN THE DECAYS 2*+p+7° AND A+p+7* 


E. F. Beall, Bruce Cork, tf D. Keefe, P. G. Murphy,! and W. A. Wenzel 
Lawrence Radiation Laboratory, University of California, Berkeley, California 
(Received September 1, 1961) 


In recent years, several authors have proposed 
theories that predict, either on the basis of vari- 
ous forms of global symmetry and the |AT| =2 
rule,*“® or on the basis of extended chirality in- 
variance,” that the asymmetry parameters in the 
decays 5+ +p+7°%(a,) and A+p+m~(a,q) should 
obey the relations ag=-a,. Here we have 


a = 2Re(S*P)/(1S1? + IP 1?), 


where S and P are the amplitudes for the two pos- 
sible angular-momentum channels in each decay. 
Other theories predict the same sign for a, and 
a,.*°*° Asymmetry measurements have shown 

la,l and la,l to be large.” Two published 
measurements of the sign of a, are in disagree- 
ment.’*44 The experiment reported in this Letter 
was designed to establish the signs and magnitudes 
of both a, and a, by measuring the polarization 

of the decay proton from £*+p+m° and A+p+m~ 
with a carbon-plate spark chamber. 






Figure 1 shows the apparatus used in the experi- 
ment. Positive pions of 1.19-Bev/c momentum 
from the Bevatron were incident upon a liquid- 
hydrogen target, producing the reactions 7* +p => 
=*+K*, £*+p+n°. During approximately one- 
third of the run, the hydrogen target was replaced 
by a block of lithium deuteride. In this case, 7* 
mesons of 1.02-Bev/c momentum produced the 
reactions 7++n+A+K*t, A+p+m~. The produc- 
tion of a + or A hyperon was indicated as in ear- 
lier experiments”® by the identification of a K* 
with a counter telescope, including detection of 
the decay of the K* in the large water Cerenkov 
counter Cx. The hollow-plate spark chamber in 
the K* telescope and the carbon-plate “proton” 
spark chamber were triggered by a coincidence 
between the K* telescope signal, the signal from 
the “proton” counter telescope that detected parti- 
cles with v/c<0.75 entering the carbon-plate 
chamber, and the pulse from a gas Cerenkov 
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FIG, 1. 


Plan view of the apparatus used in the experi- 
ment. Hyperons were produced either in the H, target 
or in the LiD target. Scintillation counters S;, S3, and 
Ss detected the K* meson; S», Sy, and Sg are scintilla- 
tion counters that detected the decay proton. The 2- 

in. water Cerenkov counters C1, Cy, and C3 were used 
to reject particles with v/c>0.75. The water Cerenkov 
counter Cx stopped the K* and detected its decay parti- 
cle; scintillation counter S4 rejected events not originat- 
ing in the target; scintillation counters U and D detected 
pions from A decay. The K spark chamber consists of 
five 1-in. hollow plates with 0.003-in. aluminum sur- 
faces. The proton spark chamber consists of four 1-in. 
hollow plates and 49 0.5-in. graphite plates sprayed 


with silver paint. Both chambers were filled with argon. 


counter which was located in the incident beam 
to detect the 7* meson. Two 90-deg stereo views 
of each chamber were then photographed. The U 
and D counters detected those 7~ mesons from A 
decay that went approximately up or down with re- 
spect to the production plane. 

All =* events with a single proton scattering 
in the carbon plates of greater than 3-deg project- 
ed angle in either view were measured. (Only 
those events with greater than 5 deg in either 
view were used in the analysis.) The measured 
K* and proton space angles and proton range al- 
lowed us to make a once-overdetermined kine- 
matical fit to the £* production and decay. The 
scattering angle and residual range were used 
to determine the analyzing power of carbon for 
each event, with the analyzing-power data which 
have been summarized by Birge and Fowler.** 
Acceptance criteria involving goodness of kine- 
matical fit, maximum and minimum scattering 
angle, and scattering within a fiducial region 
of the chamber were imposed. Out of approx- 
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imately 25000 pictures, 358 events satisfied 
all the acceptance criteria and were used as 
data. 

If the polarization of the hyperon is nonzero, 
the helicity of the decay proton is no longer mi- 
nus a.4® However, we have 


N(0)p_(@) -N(m -0)p_ (a -@) 
V V = -a@ cos@ 
N(6) +N(a -6) * re 








where @ is the angle between the proton momen- 

tum and the normal to the =-K production plane, 

py(@) is the polarization of the proton perpendic- 
ular to the production plane, and N(@) is the num- 
ber of events. Similarly, we have*® 


_ +N(r - 8)p,,(a - ~* a 
N(6) +a - 0) ieee 





where y is the angle between the proton momen- 
tum and the direction in the =-K production plane 
which is perpendicular to the hyperon momentum. 
All of the above quantities refer to the hyperon 
rest frame. We have used the approximations 
that the proton and hyperon momenta are parallel 
in the laboratory system and that the polarization 
of the proton is the same in the hyperon rest 
frame as it is in the laboratory system. There- 
fore, the likelihood functions, 





358 
L,(@) - : — -a cos@ A. —, 


and 


358 
L,,(@) - He -a cosy A . mnt: ), 


were computed. Here, A is the carbon analyzing 








power and ¢g is the angle between the plane of 
scattering and =-K production plane. The com- 
mon logarithms of Ly(a), Ly(a@), and L(q@) 
=Ly(a)Ly(q@) are plotted in Fig. 2. It should be 
noted that Ly(a@) and Ly(a@) independently imply 
that a, is positive. The combined data give a, 
=+0.75+0.17. Corrections for systematic effects 
have not been included. This value for a, agrees 
well with the lower limit for |a@,| obtained from 
a measurement of agPs by Cork et al.?° 
Because of the Fermi momentum of the target 
neutron, we could not use precise kinematical 
fitting on the A events. Therefore, the U andD 
counters were used to determine whether the 7" 
meson went up or down with respect to the K+-1* 
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FIG, 2. Common logarithms of the likelihood func- 
tions L y(a) =[];(1 -@ A; cos0; cos¢s;), L ya) =[]j(1 -0 4; 
xcosdj sings ;), and L (a)=Ly(@)L g(@) for 2*— p+ n°, 


plane, and a measurement of the polarization of 
the proton perpendicular to this plane correlated 
with the presence of a U or D count allowed a 
determination of the sign of a from the relation 


[orn “Ny boyy 
Ny +N, 





|- -a(cosé), 


with averages taken over all proton directions. 
From measurements of aap, by Crawford et 
al.,"" it is known that la, | >0.75+0.17. If it is 
assumed that |a,| >0.6, a maximum -likelihood 
analysis of our data indicates that a A is negative 
with a confidence of at least 25:1. This result 
disagrees with the early experiment of Boldt et 
al.,"* but agrees with the recent experiments by 
Birge and Fowler’® and Leitner et al.’” Thus, 
the evidence that a a is large and negative is 
how quite strong. 

The above results imply: (a) The predictions 
by Nakamura and Konuma’ that a, <0 and that 


aa >0O are in disagreement with experiment. 

(b) The phenomenological theory by Bludman,® 
which predicts a, =+aQ, is in disagreement 
with experiment. (c) The S-wave bound-pion 
model of Barshay and Schwartz® (which has been 
quoted as evidence for odd D-A relative parity”) 
is also in disagreement with experiment. (d) The 
evidence that a(=-_. a 477) is positive’® is 
strengthened. (e) A fundamental prediction of 
the global (or doublet) symmetry theories,’ 
viz., @\=-ag, is now confirmed.” 

We are grateful to the many persons who as- 
sisted in the building of the equipment, the set- 
ting up and operation of the experiment, and the 
analysis of the data. 
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PERIPHERAL COLLISIONS AND THE (3, 3) RESONANCE 
IN THE REACTION p+p+-n+p+m* AT 970 Mev 
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(Received September 13, 1961) 


As a continuation of earlier work with a diffu- 
sion cloud chamber’ at the Birmingham University 
1-Bev proton synchrotron, proton-proton scatter - 
ing has been studied with improved accuracy and 
statistics using a 9-inch liquid hydrogen bubble 
chamber.” Protons were scattered out of the syn- 
chrotron at 4° from a carbon target, and the col- 
limated beam was estimated to have an energy 
spread of +10 Mev at the bubble chamber. 3007 
events within a 1° beam spread and a restricted 
fiducial region of the chamber have been analyzed 
so far, yielding the cross sections shown in Ta- 
ble I. 


Table I. Cross sections for p-p scattering at 970 
Mev. 








No. Cross section 
Reaction of events (mb) 
(1) p+p—ptp 1554 24.4 +1.0 
(2) —d+nt 35 0.55 +0.1 
(3) —ptn+n* 1170 18.4 +0.8 
(4) —ptpr+n° 239 3.8 +0.35 
(5) —pt+ptaten 1 eee 
(6) —ptn+nt+n° 1 ~_ 
(7) —d+nt+7° 1 ar 
Either (3) or (4) 6 eee 





“These cross sections are normalized to a total cross 
section of 47.3 +1.0 mb interpolated from the many 
counter experiments. The elastic scattering cross 
section includes a correction for scanning losses at 
very small angles. This table does not include data 
from the diffusion cloud-chamber experiment, where 
the number of uncertain events was greater. 


288 


In separating reactions (3) and (4), the accuracy 
of measurement was such that it was frequently 
necessary to assume only one pion to be created 
in the interaction and to identify proton or 7* on 
the basis of bubble density. However, since only 
one case of reaction (6) was definitely established, 
and only one example of double-charged pion pro- 
duction (5) was observed, it is unlikely that the 
events classed as reactions (3) and (4) contain 
many cases of double pion production. In only 
six cases was it impossible to distinguish between 
(3) and (4). A useful check on the identifications 
is that the center-of-mass angular distributions 
were symmetrical backwards and forwards within 
statistical errors. 

Chew and Low’ have suggested that inelastic re- 
actions in which one pion is created should be 
strongly influenced at low momentum transfer 
by poles in the S matrix lying just outside the 
physical region. The physical interpretation is 
that there is a large cross section for peripheral 
collisions of one nucleon with a virtual pion in the 
field of the other, with a small momentum trans- 
fer, A, to the latter. In the present experiment 
there are two poles, since the initial state con- 
sists of two protons, and the resulting cross sec- 
tion must be symmetrical between them. 

Bonsignori and Selleri* have pointed out that, 
in p-p scattering, peripheral interactions will be 
most important in reaction (3) where the neutron 
acts as “spectator” and the proton and m* interact 
in a T= % state. They find good agreement with 
the diffusion cloud chamber results for low mo- 
mentum transfer. Using this approach, Selleri® 
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has deduced the formula 


ao _ 4f*m 
Teo ia wR (w)o(w)[a(T) + d(T, w) +c(T, w)], 


(1) 


where f* is the renormalized pion-nucleon coup- 
ling constant 0.08, m and yu are the nucleon and 
pion masses, and p; is the laboratory momentum 
of the incident proton. T is the laboratory kinetic 
energy of the “spectator” nucleon, and w is the 
total energy of the pion and the other nucleon in 
their own rest frame. o(w) is the 7-nucleon 

cross section at a center-of-mass energy w, and 
the functions R, a, and b are defined by 


R(w) = $[w* - 2w?(m? + u?) + (m? - 7)? 7, (2) 


a(T) =A,°/(A,* +4)’, (3) 
b(T, w) =A,2/(A,? +n), (4) 

with 
A,” =2mT, (5) 

and 
A? =m? +2m(T , -T) -w?, (6) 


where T; is the laboratory energy of the incident 
proton. The poles are defined by zeros of the 
denominators of (3) and (4), and are shown as 
dashed lines in Fig. 1(a). c(7,w) is an interfer- 
ence term between the two poles and is defined 
in reference 6. 

The validity of this theory is clearly demon- 
strated by Fig. 1(a), where the Q values® of (p7*) 
pairs are plotted against the laboratory energy 
of the neutron for 1164 events of the type p+p- 
p+n+a*. It may be shown that the phase-space 
factor p for this reaction is given by 


p=[R(w) /wQaT, (7) 


being a function of Q only. In Fig. 1(a) the scale 
has been distorted to allow for this in such a way 
that points should be uniformly distributed on a 
Statistical basis. The clustering of points close 
to the phase-space boundaries clearly demon- 
strates the effects of the poles and the influence 
of the well-known (3, 3) resonance at 140 Mev, 
and it accords reasonably well with the contours 
of constant cross section calculated from (1) and 
shown in the figure. 

Figure 1(b) compares the experimental values 
of do/dT with the predictions of the statistical 
model and with Selleri’s formula. The latter is 
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FIG. 1. (a) A scatter diagram of Q values of the (pn*) 
isobar against neutron laboratory kinetic energy for the 
reaction p+p—-p+n+ n*, The locations of the poles 
are shown by the dotted lines. Contours of the cross 
section predicted by Selleri’s formula are shown in 
arbitrary units from 1 to 25. (b) The neutron labora- 
tory energy spectrum compared with Selleri’s formula 
and with phase space normalized to the total number of 
events. 


in good quantitative agreement for momentum 
transfers less than 600 Mev/c (T less than 180 
Mev) and in qualitative agreement for greater 
momentum transfers. A similar agreement was 
obtained by Selleri® with the diffusion cloud-cham- 
ber data. 

Using Eq. (1), an attempt has been made to in- 
fer the 1*-p cross section using only events with 
momentum transfer to the neutron less than 430 
Mev/c (100 Mev). Events close to both poles have 
been used by calculating the momentum transfer 
to each nucleon, and taking the one with the lower 
value as the “spectator”; this amounts to folding 
Fig. 1(a) about T vaiues midway between the poles. 
Since the theory assumes contributions only from 
the one-pion exchange pole, the agreement shown 
in Fig. 2 with the known 7*-p cross section (full 
line) is quite good, and provides some justifica- 
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FIG. 2. The cross section for nr -p scattering in- 
ferred from Selleri’s formula using only those events 
in which the momentum transfer to the neutron is less 
than 430 Mev/c. The full line shows the cross section 
determined from direct rt=p scattering experiments. 


tion for similar studies of the 7-7 cross section 
from inelastic 7-p scattering.” However, it is 


found that including events with momentum trans- 
fers from 430 Mev/c to 530 Mev/c reduces the 
agreement of Fig. 2, so that other considerations 
are necessary in this region. A guide to the mag- 
nitude of interference effects is that in the region 
of momentum transfer up to 430 Mev/c, interfer- 
ence of the amplitude due to the second pole 
through the term c(7,w) in (1) amounts to 5 to 
10%, according to the Q value; it has been al- 
lowed for in plotting Fig. 2. The original pro- 
posal of Chew and Low’ was to infer the 7-p 
cross section for a number of strips of A? and 

to extrapolate to the pole, when interference due 
to other poles disappears. This has been done 

by Smith et al.,° but there are insufficient events 
in either their or our experiment to carry out 
this procedure except in wide intervals of Q. 


The comparison of the present results with the- 
ory shows clearly that most of the interactions 
arise from peripheral collisions, thus accounting 
for the success of the isobar model of Linden- 
baum and Sternheimer’’® in p-p experiments. 
For those events with momentum transfers great- 
er than 600 Mev/c, the Q-value distribution is, 
however, in good agreement with the statistical 
model. For 7-p scattering, where peripheral 
collisions and isobar formation have different 
roles, the situation is less clear.*?*° 
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REMARK ON THE ALGEBRA OF INTERACTIONS 


A. Pais* 
CERN, Geneva, Switzerland 
(Received September 8, 1961) 


It is the purpose of this note to point out some 
remarkable similarities between the structure of 
interactions and the algebra of octonions. 

To characterize briefly this strange algebra, 
we recall some familiar properties of quaterni- 
ons. 


X =X ,e 


o°0**a*a (1) 


(summation over a =1, 2,3) is a quaternion if 


(2) 


ee +e = -26 
a 


B*°p°a ap’ “ap apy y' 


Egy is totally antisymmetric; €,,,=1. If the 
“components” X9, Xq are real, define 


(X,Y)=3(XY +¥X), X=X e -e xX , (3) 
00 aa 


N(X) =(X,X). (4) 


In this real case the norm of the product XS of 
two quaternions satisfies 


N(SX) =N(XS) =N(X)N(S), (5) 
from which by associativity 
N(S74XS) =N(X). (6) 
A null quaternion is defined by 


X=0: X,=X_ =0. (7) 
0 a 
From Eqs. (2) and (3), (€q, €g) =6qg so that the 
three ey are like orthogonal unit vectors. 
Quaternions share with real and complex num- 
bers the property (5). There exists only one oth- 
er number system’ (with two variants) for which 
Eq. (5) is true, namely, the octonions. The first 
variant satisfies in fact all Eqs. (1)-(7)—with two 
changes of meaning, however: (a) a,8,+++ =1, 
oo, 75 (b) Eqpy is totally antisymmetric and 
equals +1 for 


(aBy) =(123), (145), (167), (246), (275), (365), (374). 


(8) 


Thus one generates octonions from our starting 


quaternion by introducing a further unit vector 
e4, (€a,e4)=0, a=1,2,3. Then by vector prod- 
uct formation 2e,=e, Xé€,, 2€,=e,Xe,, -2e,=e, 
xe,. Equations (2) and (8) show that octonion mul- 
tiplication is nonassociative and hence does not 
yield a group.” However, any two octonions gen- 
erate a group; hence Eq. (6) is true as well. 

The continuous automorphisms of the €g form 
the group G,. In addition we note the following 
two discrete operations. 


Pi: C9 > 0, (9) 


P,: ey eye @=4,5,6,7; i=(-)™. (10) 
P, does not, P, does change the multiplication 
table (8). P, adjoins’ to an X a “split” x‘), x(s) 
=Xret+iXjm, where Xre, Xjm are orthogonal. 
Define N’(x‘S)) =N(Xre) -N(Xim); then 


(s) (s) 


nix )y'S)) nex 


N'Y”), 
an alternative version of Eq. (5) (second variant). 

q-octonions have q-number fields as compo- 
nents. Here too one can define an inner product 
by taking the Hermitian average on the right side 
of Eq. (3). Let X be a qg-octonion, S a real c- 
octonion. One shows that Eqs. (5) and (6) are 
again true. Thus if N(S)=1, then N(XS) =N(SX) 
=N(S~XS) =N(X), i.e., one can define octonion 
gauge transformations of qg-octonions. 

Let B=Bog -i€gBg, M=MgQ-i€qgMq, where the 
components of B(M) are spin 3(0) fields. The 
equation* 


(y8 +m)B +iGy,BM =0 (11) 


expresses, by Eqs. (2), (7), and (8), the strong 
interactions with a high symmetry provided B, 
=A, 2¥9B,=2++2", i2”B,=2°-2*, B°=2°, 
i2!B =n+=°, -2B.=p+= , i2”B,=p-= , 
-2”B,=n-=°. Put B=B(A,2,N,=). Then® 


M =B(o,7,K,K°). (12) 


The M components are eigenstates of charge con- 
jugation. The “full symmetry” is: global sym- 
metry of the type® G~ and moreover K -coupling 
strengths = minus 7-nucleon coupling. 

Thus the octonion calculus which involves sets 
of eight “equivalent particles” automatically pro- 
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duces all needed selection rules plus those un- 
wanted extra ones implied by too strong a sym- 
metry. We now observe that we can adjoin to B 
two other octonions, namely B“ =P,B and B® 
=P,RB, where R=(e,-0,4, €4~ -€7, C5 > €5, Cg > 
-e,) and is in G,. Replace in Eq. (11) BM by 

[B +c,B+c,B™]M and the following happens: 

In the right qualitative way, the baryon masses 
split as 8=1+3+2+2, the meson masses as 8 
=1+3+4. Each of the three couplings separately 


have equivalent full symmetry, denoted by F, F™, 


F™, respectively. Some partial symmetries are: 
F +F® =doublet approximation, F +F ” =Beh- 
rends-Sirlin scheme.” Note that three is the 
minimal number of clashing full symmetries 


which breaks the degeneracies.*® The usual charge 


operator 
Q=T,+Y, (13) 


is the sum of the third components of isotopic 
spin and hyperspin and has the following curious 
property. With respect to each of the F’s sep- 
arately, Q itself is isomorphic to a third compo- 
nent of angular momentum. (This is not true for 
any other nontrivial linear combination of T, and 
Y..) 

The above example of the lifting of degenera- 
cies is not unique. In particular, the present 
calculus is not tied to the strange particle pari- 
ties.° 

Consider next the transformations 


B-S™"*BS, M+S™MS, (14) 


with S a real c-octonion, N(S)=1. By application 
of G,,S can be brought into the canonical form S 
=expe,. For infinitesimal é, the transformation 
(14) is not in G,. Even so, Eq. (11) with G=0 is 
invariant under Eq. (13) as S, B generate a group. 
For G#0, Eq. (11) becomes (y8+m)B+S[(S“BS) 
x(S“MS)]S“*=0. For “infinitesimal” = one finds 
that now BM ~ BM +“weak” interaction with AT 


=3. Thus e, acts as the spurion. For the present 
we do not discuss the parity structure of weak non- 


leptonic interactions so generated. 

If the octonion algebra envelops in some sense 
the structure of interactions, one may ask how 
the leptons could enter. Here the electric charge 
should form one bridge. Especially if two neu- 
trinos v,,v, exist, it is interesting to contemplate 
the possibility that, rather like M, the lepton is 
a self-charge conjugate octonion L. L can be so 
constructed that, where (B, Qy,B) is the electric 
baryon current, so (L, Qy,L) is the electric lep- 
ton current. In treating neutral leptons consider- 
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able arbitrariness exists, however, partly con- 
nected with the group SU(3) (subgroup of G,) of 
automorphisms which keep one é, fixed. Lep- 
ton conservation induces insufficient restrictions, 
Some partial results follow. 

Using doublet language,*° write B as B(N,,N,, 
N,,N,). Put’? 


L(w) =2-B(wn,°, -wn,°, ng, 4), (15) 


no(@). m= (3) 


w is a 2X2 unitarity matrix; (wtr,w)? = -T,. As 
for B, one checks that Q is indeed the L electric 
charge operator. 

Define current octonions” J=LyL, j=Byl +ay,) 
xB. The inner product, 


(J, j) * Soda 
describes leptonic transitions, as follows. a=0,3: 
AS =0, neutral; a=1,2: AS=0, charged; a =4,7: 
|AS| =1, neutral; a=5,6: |4S| =1, charged. One 
finds J,=J,=J,=0 for any w as defined. Hence 
there are no neutral |AS| =1 transitions. Let J” 
correspond to w=itT,. One finds J,” =J,%=0. 
Thus (J +j, J” +7) describes the universal AS 
=0 interaction.’* Let J refer to (a) interchange 
of v, with v;°, v, with v7°, (i#j) =1, 2; (b) w=ir,. 
One finds J,” =J,=0, while now J,”, J, are 
lepton-conserving |AS| =1 charged currents. The 
two choices for (i, j) are related to the neutrino- 
flip question.** It is perhaps a good aspect that the 
absence of |AS]/ =1 neutral leptonic follows from 
a specific algebra. It is perhaps a bad aspect 
that the synthesis of AS =0 leptonic and |AS|=1 
is so far not unique. (The use of different L’s 
is not unlike the use of different B’s for strong 
interaction asymmetries.) 

It should be noted that all results stated here 
could have been written without using octonions 
at all. The formal structure here described 
seemed sufficiently intriguing to communicate, 
however. Yet the close connection between oc- 
tonion algebra and some aspects of the interac- 
tions may be nothing but a strange coincidence. 
To be more than that, this algebra should play 
a dynamical role.*® 

The author is much indebted to CERN for its 
hospitality and financial support. He also wishes 
to acknowledge a consultation with I. Ching. 





*Permanent address: Institute for Advanced Study, 
Princeton, New Jersey. 
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‘Hurwitz’ theorem. For early developments see L. 
Dickson, Ann. Math. 20, 155 (1918). For both variants 
the components are real (real octonions). 

*But rather a so-called quasi-group where the associ- 
ative law is replaced by a certain weaker law. For the 
present it suffices to state that for octonions this law is: 
The associator X(YZ)-(XY)Z is totally antisymmetric 
inX,Y,Z. 

3In this case and also for the g-octonions mentioned 
below, we have no longer a division algebra. This does 
not matter as we need not take inverses of wave fields. 

‘For the purpose of exposition we consider only tri- 
linear interactions and leave aside the question whether 
strong interactions are generated by gauge fields. 

59 is a T=S =0 meson conjectured by many authors. 
0=0 is not necessarily excluded. K=(K*,K°); K© 
=(-K°,K~). 


8In a terminology employed elsewhere; A. Pais, Phys. 


Rev. 122, 317 (1961). 

tk, Behrends and A. Sirlin, Phys. Rev. 121, 324 
(1961). 

‘Example: B® =R’B; R’ =(eg > -e€q), @ =4,5,6,7 and 
is inG). Even with all these couplings there remains 
a simple 7-coupling constant relation: (NN7)=-(Z=7) 
=(Z=7). Note: P, corresponds essentially to BM— MB. 

‘Example: Write the above B as B(N,A,=Z,=), define 
B'=B(N ,cy¥sA,2 ,Co¥s=), and replace in Eq. (11) BM by 
B'M. This spreads the masses due to space parity ef- 
fects. The curious property of Q is now lost. The ex- 


ample in the text may be said to generate splits due to 
two “isotopic parities.” 

‘Tn the notation of reference 6. As is done in part of 
the quoted paper, the doublets are used as a mathemat- 
ical device without necessarily insisting on any mass 
degeneracy. 

‘Ic — charge conjugate. v;=3(1 +Y5)U;- ¥; is a mass- 
less spinor. It is possible but unattractive to put 1; 
=v). T =transpose. 

The four-vector index is suppressed. J will always 
have the correct V-A combination due to the very def- 
inition of »;. The jg, @=4,5,6,7, have AT=4, Also 
AT =4 currents can be constructed, namely by using 
G» generators; see reference 7. 

Including a neutral 4S =0 current J3°" without ve 
terms. To such acurrent there is so far no objection; 
see S. Bludman, Phys. Rev. 115, 468 (1959). (In the 
present case the Wu and €e terms moreover conserve 
parity.) For a first attempt to tie lepton phenomena to 
a group, see S. Bludman, Nuovo cimento 9, 433 (1958). 

Kg. Feinberg, F. Giirsey, and A. Pais, Phys. Rev. 
Letters 7, 208(1961). Another ambiguity exists due to the 
fact that j need not necessarily be the same for AS = 0 
as for AS=1. 

SNote the possibility of linearized octonion wave equa- 
tions. Let X be an octonion field, O=eypy+iepmp+mo, 
w=1-4, p=5-7, mp mass parameters. Put OX=0. Then 
0(OX) =(00)X =0. Hence we get a standard wave equa- 
tion for mass (m,° + mo’), for each component of X, 
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